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Abstract 

A combination of spectral and moments analysis of the continuous 
X-ray flux data is used to show consistency of statistical properties 
of the powerful solar flares with 2D BTW prototype model of self- 
organized criticality. 
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1 Introduction 



The solar flares are the manifestation of a sudden, intense and spatially con- 
centrated release of energy in the solar corona. According to the Parker's idea 
PP stochastic photospheric fluid motions shuffle the footpoints of magnetic 
coronal loops. The high electrical conductivity of the coronal gas implies 
that the magnetic field is frozen-in, so that the subsequent dynamical relax- 
ation within the loop results in a complex, tangled magnetic field, essentially 
force-free everywhere except in numerous small electrical sheets which form 
spontaneously in highly-stressed regions. As the current within these sheets 
is driven beyond some threshold, reconnection sets in and releases magnetic 
energy, leading to localized heating. In this picture there exists a separa- 
tion of timescales between energy input to the system (minutes to hours, 
for photospheric motions), and energy release (seconds to minutes, for re- 
connection and thermalization under coronal conditions). It means that the 
coronal loop is slowly driven by footpoint motions. Thus, the coronal field 
does meet the requirements for the appearance of self-organized criticality 
(SOC): self-stabilizing local threshold instability; open boundaries; and sep- 
aration of time scales between driving and avalanching. Following to this 
picture the authors of Ref. [2] suggested that SOC could describe the main 
features of the hard X-ray flares with avalanches of the reconnection events as 
the flares (see for further development and references therein). Power 

law probability density functions (PDF) observed for peaks, waiting and du- 
ration times of the X-ray flares present the main experimental support of 
this hypothesis. This experimental information, however, turned out to be 
non-sufficient to distinguish between the SOC and fluid (MHD) turbulence 
processes mixed in the solar corona jH]- There exist substantial differences 
between SOC and turbulence mechanisms for explanation of the power law 
PDF's jHj- Such differences can be ascribed to the lack of long time corre- 
lations for the avalanches in the monoscale SOC models. The distinction 
between the two mechanisms is however lost in the 2D BTW (Bak-Tang- 
Wiesenfeld) model |7],jH] which obeys a peculiar form of multiscaling for the 
PDF's of several avalanche measures 0-^1] . The 2D BTW model has long 
time correlated bursts and other turbulence-like intermittent properties if 
studied at the wave time scales. It is suggested in [TT] that such SOC models 
would be candidates to describe the solar flares. It should be also noted, that 
recent investigations P2]>E3 snow that the lack of exponential distribution 
for the waiting times between the solar flares (which is mentioned in Ref. jB]) 
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cannot be considered as a reason for discarding the SOC models as under- 
lying dynamics. The long-range correlations of different nature can result in 
the non-Poissonian behavior of the waiting times for the SOC systems. The 
long-range correlations for the solar flares dynamics can have both internal 
and external reasons (a correlated drive related, for instance, to planetary 
motions affecting the Sun dynamics [T5]). The problem now is: How can 
one distinguish between fluid (MHD) turbulence and the SOC models with 
long-time correlations in the observed flares? We consider continuous X-ray 
solar flux measured during a remarkable two-week period from 9 July to 23 
July 2000 year. This period can be called as a flares festival: 98 C-class, 39 
M-class and 3 X-class flares. For example, in analogous period of 1997 year 
the only one low-power B-class flare was observed. In the minute time scale 
we have more than 21000 data points presumably dominated by the powerful 
X-ray flares. We will use a combined spectral-moments analysis of the signal 
to extract information about underlying physical mechanisms. 



2 Data and analysis 

The data was obtained by GOES-10 satellite of the National Oceanic and 
Atmospheric Administration USA (NOAA) providing the kind of continuous 
monitoring of the whole-sun X-ray flux, J, for the 0.5-to-4 A wavelenght band 
(1-minute averages). The data are given in Watts /cm 2 sec. 

We believe that a mix of different complex mechanisms produce the to- 
tal X-ray flux. Therefore, characteristics of the raw signal generally do not 
exhibit some recognizable properties, which could be used for identification 
of the underlying physical processes. One should use some "threshold" tech- 
nic to obtain the recognizable features. On the other hand, even if such 
features will appear, then a single certain characteristic, such as spectrum 
or PDF, cannot univocally indicate a dominating mechanism. For instance, 
different mechanisms (including SOC and fluid turbulence) can produce the 
same power-law spectrum To resolve this problem we will analyze two 
different physically meaning fields: the flux, /, itself and its so-called local 
dissipation rate [TT].[r3]: 

e t = J2(h + i-Ik) 2 /t (1) 
fc=i 
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Statistical properties of such gradient measure are generally independent 
on statistical properties of the original time series This fact allows us 
considering simultaneously the two statistically independent fields: I and 
e = (A/) 2 , produce a reliable identification of the underlying physical mech- 
anisms. The name of the e measure comes from fluid turbulence, where the 
measure indeed is directly related to the dissipation processes |T5] . More- 
over, in our case the e field seems to be also directly related to dissipation 
(see Introduction). Therefore, we will perform the threshold discrimination 
of the data set using just the dissipation field = (h+i — h) 2 - Then, we 
return to the /-field to look at spectrum of / after excluding the data points 
with the large (above the threshold) dissipation. The recognizable power-law 
spectrum of / appears at low-frequencies at the threshold value of the dissi- 
pation field e ~ 1CT 12 W 2 /cm 4 s 2 . The spectrum is shown in figure 2 (figure 
1 shows the low-frequency spectrum of the raw flux / for comparison). 

The straight line is drawn in figure 2 in order to indicate the power-law 
spectrum 

E{oj)~u)- 1 (2) 

in the log-log scales. For the discriminated data set (with the spectrum (2) 
at the low frequencies) we calculated moments of the local dissipation rate 
(1): (ef). For the multiscale SOC (in particular for the 2D BTW model [TT] ) 
appearance of the multiscaling 

(e?) ~ (3) 

is expected. Figure 3 shows the moments against t (in log-log scales) and the 
straight lines are drawn in the figure in order to indicate the multiscaling (3). 
Slopes of these straight lines (best fit) provide us with multiscaling exponents 
fi p , which are shown in figure 4 as circles. The time interval where we observe 
the multiscaling of the dissipation field is: 10 min <Tt < 1000 min. 

3 Discussion 

The BTW prototype model of self-organized criticality is defined on a square 
lattice (L x L) box jHJ. On site i Zi = 0, 1, 2, . . . is the number of "grains". 
If z% < 4, Vz, the configuration is stable. Grain addition to a stable con- 
figuration is made by selecting at random a site k where z^ — > Z}~ + 1. If 
then Zk > 4, toppling occurs, i.e. z^ — > Zk — 4, while each nearest neighbor, 



4 



/, of site k gets one grain (zi — ► zi + 1). If is at the border grains are 
dissipated. The toppling of site k may cause instabilities in the neighbors, 
leading to further topplings at the next microscopic time step, and so on. 
Thus, an avalanche made by a total number s > of topplings occurs before 
a new stable configuration is reached and a new grain is added. After many 
additions the system reaches a stationary critical state in which avalanche 
properties are sampled. 

A key notion in the approach to the BTW is that of wave decomposition 
of avalanches JH] • The first wave is obtained as the set of all topplings which 
can take place as long as the site of addition is prevented from a possible 
second toppling. The second wave is constituted by the topplings occurring 
after the second toppling of the addition site takes place and before a third 
one is allowed, and so on. The total number of topplings in an avalanche is 
the sum of those of all its waves. The wave size has a pdf satisfying finite 
size scaling P w (s,L) ~ s~ Tw f w (s/L Dw ), where r w = 1, D w = 2, and f w is a 
suitable scaling function. 

The multiscaling of the dissipation field 

t 

k=l 

was observed in the Ref. ^1] for the 2D BTW model and the multiscaling 
exponents for the BTW dissipation calculated in JT] are given in figure 4 as 
triangles. One can see good correspondence between the exponents calculated 
in jTTj and the observed exponents. 

The dissipation threshold e ~ 10~ 12 W 2 /cm i s 2 1 for which the solar X- 
ray data agree with the 2D BTW model both for the flux spectrum and for 
the dissipation multiscaling, corresponds to transition from the low intensity 
M-class flares to the high intensity M-class flares. At present time we do 
not know why just at this level of the dissipation rate the SOC (presumably 
BTW) mechanism overcomes all other mechanisms and appears so clear in 
the X-ray flux data. In any way, we can see that the multiscale SOC takes 
place there and even appears as a dominating mechanism at a certain stage. 

The authors are grateful to the NOAA for providing the data. 
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Figure 1: Spectrum of the X-ray flux raw data I (in log- log scales). 
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Figure 2: Spectrum of the X-ray flux data set, /, after exclusion of the points with the 
local dissipation rate e > 1CT 12 W 2 /cm 4: s 2 . The straight line is drawn to indicate the 
power law (2). 
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Figure 3: The local dissipation rate moments (p = 2, 3, 4) against t in log-log scales. The 
data set is the same that was used to obtain figure 2. The straight lines (the best fit) are 
drawn to indicate the multiscaling (3). 
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Figure 4: The multiscaling exponents fi p (3) extracted as slopes of the straight lines in 
figure 3 (circles) . Triangles correspond to the multiscaling exponents calculated in ^I] for 
the 2D BTW model. 
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